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Abstract 


The theory of aging of a poorly soluble precipitate such as AgBr by dissolution- 
reprecipitation according to the higher solubility of small particles compared to 
that of large particles is developed quantitatively. Equations are given for the 
temporal change in the size distribution function, the average particle size and the 
particle number for diffusion-controlled and reaction-controlled mass transfer. 


1 Problem statement 


According to the GIBBS-THOMSON formula, small droplets have a higher vapor 
pressure or a higher solubility in a surrounding liquid phase of a second substance 
than large drops. In a multidisperse emulsion, matter is transported from the 
surface of small drops to the surface of large drops, i.e. large drops grow at 
the expense of small drops. Analogous processes can be found in multidisperse 
suspensions of crystalline substances. As an example, the grain enlargement of 
AgBr precipitation that occurs in the production of photographic emulsions is 
mentioned, which is known in the literature as OSTWALD ripening. 

The task of this work is to develop the equations that describe these processes 
quantitatively. This assumes that agglomeration of particles does not occur or 
is only a side effect of minor importance. The initial state is characterized by a 
certain distribution function of the size of the particles. Therefore one can ask, how 


1. the distribution function, 
2. the average particle size and 
3. the total number of particles 


changes over time. The same task is handled, among others, by TODES, [15], 
TODES and KHRUSHCHEY, {14] as well as by LirsHiTz and SLYOZoV, {9, 10]. The 
results of these authors are referred to below. 


2 List of Symbols 


Co Saturation concentration for the limit case of 
infinitely large particles 

Ce Saturation concentration for particles of radius r 

é. Concentration at the surface of particles of radius 
r 

Cc average concentration in a system of particles of 
different sizes 

F Diffusion constant 

f (r,t) Distribution function of particle number as a 


function of radius r at time t according to the 
definition in Equation 13 

gp(t), hp(@) Auxiliary functions for diffusion-controlled (D) 
and reaction-controlled (R) mass transfer, cf. 
Equation 53 and Equation 54 

gr(t),he(e) Auxiliary functions for diffusion-controlled (D) 
and reaction-controlled (R) mass transfer, cf. 
Equation 73 and Equation 74 

k; Constant of the phase boundary reaction, cf. 
Equation 19 


List of Symbols 


YD; YR 


Number of moles in a single particle 

change over time in the number of moles in a 
single particle 

effective particle radius, cf. Equation 3 
change in particle radius over time 


Radius of a particle in equilibrium with the 
surrounding solution 


mean radius in a heterodisperse dispersion at 
time t = 0 

Gas constant 

Temperature 

Time 

Molar volume 

Total number of particles per unit volume 


Total number of particles per unit volume at 
time t = 0 


Auxiliary parameters for diffusion-controlled (D) 
and reaction-controlled (R) mass transfer, cf. 
Equation 58 and Equation 78 

Scattering of the particle size distribution at time 
t= 

Scattering of the particle size distribution at time 
t 

-z, dimensionless measure of particle size 
Interfacial tension 

stoichiometric factor 

Time constants of dissolution-reprecipitation for 
diffusion-controlled mass transfer, cf. 

Equation 46, Equation 60 and Equation 67 
Time constants of dissolution-reprecipitation for 
reaction-controlled mass transfer, cf. 

Equation 48, Equation 80 and Equation 88 
Volume of the disperse phase per unit volume of 
the dispersion 


3 Thermodynamic equations 


3 Thermodynamic equations 


For the vapor pressure p of a droplet of radius r we have the well-known GIBBS- 
THOMSON equation: 


Pr = Po° erat ) (1) 
where po is the vapor pressure for the limiting case of a flat liquid surface (r = oo), 
a is the surface free enthalpy, V,, is the molar volume of the condensed phase, 
& is the gas constant and T is absolute temperature. An analogous formula 
applies to the solubility of small droplets and small crystals in a surrounding liquid 
phase. In order to obtain a formula that can also be used for crystals, [2, 13], it 
must be assumed that the crystal which is present is in the so-called equilibrium 
configuration according to GIBBS-WULFF, corresponding to a minimum of the 
total surface energy at a given mass or number of moles of the crystal. Then the 
exponential function 

OV, OO 

1er Z a 
arises, where o is the mean free surface energy corresponding to the equilibrium 
configuration of a particle and ge is the quotient of surface change and volume 
change when a particle increases in size while maintaining the equilibrium configu- 
ration, the presence of which is assumed in all following considerations. In case of 
a spherical particle: 


4 
(v = gtr O= tar”) 
the following holds: 
wo () 


which results in Equation 1. 
Introducing the concept of the effective radius of a particle according to the 


definition: av 
=a ones 
r 50 (3) 


one can therefore write for a non-dissociating substance instead of Equation 1: 


Cr — Co . ent 5 (4) 


where Cg is the limit of the solubility of an infinitely large particle and c, is the 
solubility of a particle with effective radius r = 2 - ( or). In general, the exponent 
is small compared to 1 and the exponential function can therefore be developed in 
form of a TAYLOR series. Thus it follows: 


cr ey: |i + (5) 


3 Thermodynamic equations 


This approximation fails for very small particles for which oo > 0.3, i.e. for 
particles whose radius is of the order of 3 x 10~’cm or smaller. However, the 
behavior of such small particles is generally practically irrelevant to the overall 
course of phenomena. Therefore, the approximate equation according to Equation 5 
is used throughout the following. 

In the case of dissociating substances (especially salts) of the general formula A,Bg, 
an equation corresponding to Equation 4 applies for the solubility product Lp: 


20Vm 


Lp, = Lp): e727. (6) 


In the case of complete dissociation and the presence of equivalent amounts of A 
and B in the solution, one can calculate the gross solubility c (moles of A,Bg per 
unit volume) by the individual concentrations ca and cpg of A and B respectively 
as follows: Se 2% 
B 
Cc SS a ee ee 7 
A= (7 
Using the defining equation of the solubility product 


Lp. =c%c8 (8) 


it follows from Equation 6 and Equation 7 at total dissociation 


20Vm 


cOt8 = Ot eraP (9) 


Using the series expansion for the exponential function we get from Equation 9: 


20 Vm ; 
io 5 1 | ; f —_ = —, 1 
nem (a+ 8)-rAT : a B 


In many cases, component A or B is present in a relatively large excess in the 
solution. If there is an excess of B, cg is given as a practically constant value. Then 
one usually defines solubility as the total concentration of the minority component 
A in the solution at equilibrium and obtains from Equation 6 and Equation 8 by 
series expansion of the exponential function: 


C= Cy f+ zd | (11) 


An analogous formula for the solubility with $ instead of a is obtained in the case 
of an excess of A measured as equilibrium concentration of B. 
In general one can write: 


20Vin 
vr kT \ 


c= eo: |l+ 


4 General equations for the time change of the distribution function in a 
polydisperse system 


In Equation 12 vy = 1 for undissociated substances. Even with dissociating 
substances, v = 1 if there is an excess of B and the formula ABg,/, or if there is 
an excess of A and the formula Ag /gB is used as a basis for calculating the molar 
volume V,,,. The above notation of the formula is assumed in the following, as this 
is also the easiest way to formulate the diffusion and reaction equations. Only in 
the case of equivalence of A and B in the solution should v = a+ {§ be set based 
on the formula A,Bg to calculate the molar volume V,,. 


4 General equations for the time change of the 
distribution function in a polydisperse system 


To describe a multidisperse system of particles distributed in a second phase, the 
distribution function f(r,t) is introduced, the defining equation of which is: 


Gf) ie MN (rir + Ar, t) 


Ar=0 Ar : ( i) 


where .V(r,r + Ar,t) is the number of particles per unit volume of the disperse 
system with radii between r and r+ Ar at time t. 

The number of particles per unit of volume that reach and exceed a certain radius 
through growth per unit of time is equal to the product of the distribution function 
and the change in radius of the individual particle over time, i.e. equal to f - a 
The change in the distribution function over time for a given value r then results 
as the negative value of the divergence! of f - o This results in the differential 
equation already given by LIALIKOV, PISCOUNOVA, CHIPILOV and CERDYCEV, 


[8]: 
af acfi) 
Ot Or - 
For the total number Z of particles per unit volume of a dispersion it follows from 
Equation 13: 


(14) 


Le [ fitydr. (15) 


The time decrease in Z is given by the number of particles that approach the radius 
r = 0 per unit volume of the dispersion and per unit time, i.e. by the limit of 
—f-r for r=0. Thus it follows: 


-F=-5(f tenar) = =r. (16) 


'The divergence is an operator, which takes in the vector-valued function, i.e. the distribution 
function over time defining this vector field, and outputs a scalar-valued function 2" measuring 


dt 
the change in radius of the particle at each point. 


5 General equations for diffusion-controlled and reaction-controlled mass transfer 


An evaluation of Equation 14 and Equation 16 requires knowledge of the change in 
particle radius r over time. The corresponding equations can be found in section 5. 
Another important general equation is obtained by calculating the fraction y of 
the disperse phase, i.e. the volume of the dispersed particles per unit volume of 
the dispersion. Using the distribution function this yields: 


oO 4 3 
y= ~ . f (r,t) dr. (17) 


The temporal change in the dissolved fraction corresponding to the temporal increase 
in the mean particle radius is generally negligible. As long as the dispersed substance 
is neither formed nor consumed by chemical reaction, y is time-independent and 
so is the integral on the right-hand side of Equation 17. This relationship will be 
used several times below. 


5 General equations for diffusion-controlled and 
reaction-controlled mass transfer 


The time change 7 in the number of moles of a spherical particle of radius r 
of a non-dissociated substance is equal to the negative value of the product of 
the diffusion coefficient Y, surface area 4rr? and the concentration gradient on 
the surface. If the concentration in the environment is significantly smaller than 
the concentration within the particle, then the concentration gradient for quasi- 
stationary conditions is one where Cc’. is the concentration at the particle surface 
and c means the average concentration of the solution, [11, 16]. So one obtains: 


Le eee ee eee (18) 
dt 

This assumes that the distance between the next neighboring particles is large com- 
pared to their radius. With regard to Equation 17, this requirement is equivalent 
to the condition y < 1 corresponding to a sufficiently small volume fraction of the 
dispersed phase. 

The above equation also applies to dissociated salts AjBg if, in the case of equiva- 
lence, the solubility is defined according to Equation 7. Furthermore, Equation 18 
also applies to an excess of B, if the calculation of the number of moles is based on 
the formula ABg;, and with ¢ there is understood the concentration of A in the 
solution and ZY means the diffusion coefficient of A. If there is an excess of A, the 
formula Aj/gB should be used. 

The mass transfer at an interface between a uniform liquid phase and the gas phase 
or between two liquid phases without the presence of additional surface-active 
substances is generally hardly inhibited, so that local equilibrium at the phase 
boundary can be expected. As a limiting case for the latter one can set that ci. = c,, 


5 General equations for diffusion-controlled and reaction-controlled mass transfer 


where c, is given by Equation 12. 

The transfer of ions from a crystal into the surrounding solution and also the 
reverse process of the transfer of ions from a solution into a crystal can also be 
practically uninhibited, e.g. in the case of AgCl or KC1Ou, [4, 3]. In other cases, 
e.g. BaSO,, noticeable deviations from equilibrium can be observed, since the 
hydration or dehydration of higher-value ions in particular can have a significant 
activation energy, [4]. In order to describe the processes also in such systems, it 
is assumed that the rate of mass transfer at the phase boundary per unit area is 
equal to the product of a rate constant and the difference between the equilibrium 
concentration c, and the concentration c’. at the surface. For the change in the 
number of moles of a particle with surface area 47r? one gets: 


n= —4ar*k -(c. —¢.), (19) 


where k means the rate constant of mass transfer. 

The relationships for nearly spherical, e.g. cube-shaped crystals are analogous. 
However, if the effective particle radius introduced in Equation 3 is used, additional 
numerical factors of the order of magnitude 1 appear in Equation 18 and Equation 19. 
This will be ignored below. The results of the following calculations therefore 
only have semi-quantitative significance for non-spherical particles. The special 
approaches for the growth of long needles or thin plates, [16], can be ignored 
here, since according to section 3 crystals are assumed to be in the equilibrium 
configuration. 

By equating the right-hand sides of Equation 18 and Equation 19 one obtains an 
equation for calculating c).. By solving for ci. we get: 


kreop + Be 
(= ee 20 
a kr+Q ce 
Introducing Equation 20 in Equation 18 or Equation 19 yields: 
Anr?k D 
ts 2 eee ee 21 
a kr+@ ace) ey) 
For the time change in the volume ar of a single particle we get: 
d (4rr? 
a(S) <tr tVn (22) 
where V,, means the molar volume of the dispersed substance. 
From Equation 21 and Equation 22 it follows that: 
tiVs be D Ve 
— = — . po “ 23 
Gare kr+Q@ ee) oe 


5 General equations for diffusion-controlled and reaction-controlled mass transfer 


The following limiting cases are to be noted: 


1. If kr >> JY, then Cc. = c, corresponds to practical reach of equilibrium on the 
drop surface, i.e. the mass transfer is controlled by diffusion. It therefore 
follows from Equation 23: 


pt NE eri, ARS, (24) 


r 


2. Conversely, if kr < Y, then c. = c corresponds to the practically reached 
diffusion equilibrium between the surface and the interior of the dispersing 
liquid, i.e. the mass transfer is reaction-controlled. It therefore follows from 
Equation 23: 

r=-kV,-(¢--—c), i kr<@. (25) 


The change in the concentration of solute over time is equal to the negative value of 
the sum of the changes in the number of moles over time of the individual particles 
per unit volume of the dispersion: 


anf He) Fetar, (26) 


If the amount of substance contained in the dispersed droplets or crystals is 
significantly larger than the amount of substance dissolved in the dispersing liquid, 
it can be assumed that the amount of substance that goes into solution from the 
small particles per unit of time is practically equal to the amount of substance 


deposited on the large particles and in comparison to this, under quasi-stationary 
conditions, the change in the amount of substance oe in the solution corresponding 
to an increase in the average particle radius is vanishingly small. Substituting 
Equation 12 and Equation 21 into Equation 26 for Ge = 0 gives a conditional 


equation for the concentration c inside the dispersing liquid with the result: 


ca le fle Bs on 
co VAT fo (r,t): Bas 
Inserting Equation 12 and Equation 27 into Equation 23 gives: 
7 _ 20eVin kQ i i So F(r,t)- ed (28) 
vrkTkr+Q fed) 


From Equation 28 one obtains for the radius r* of a particle that is in equilibrium 
with the ambient concentration and for which accordingly r = 0: 


i) r2 dr 
k to f(r, t) ; ao 


ie f(r,t)- 3 


(29) 


6 Equations for systems in which the initial distribution function is given by a 
narrow Gaussian curve 


Using Equation 29 one can write Equation 28 in the simplified form 


. 20eV2 kG r 
= vm : —1)}. 
vrk&Tkr+Q Ga ) ey) 
Inserting Equation 30 in Equation 14 gives: 
Of(r,t) 2aaVZ, Ol kD a flr,t) 
atCO* AT Or meg ee r et) 


The following limiting cases are to be noted: 


1. If kr >> ZY (with 7 as the mean droplet radius), then for the majority of the 
droplets one has diffusion-controlled mass transfer. In this case Equation 29 
and Equation 31 are simplified. r* becomes equal to the mean radius 7: 


rae fe Sher 
i Ts, t) dr 
For kr > Y it further follows from Equation 31: 


Of (r,t) 20 Doon 0 one) . 


r2 


Ot vr RT ar 


r 


= (32) 


(33) 


2. Ifkr < QY, then for the majority of droplets one has reaction-controlled mass 
transfer. For this, Equation 29 and Equation 31 simplify to: 


— des fir) ar de 


= [o°e) ) 34 
i fUyterdse a 

OFRO)..- 2akeV2 O [(r—r*)- f(r,t) 
at ert RT Or iG | , oy 


6 Equations for systems in which the initial 
distribution function is given by a narrow Gaussian 
curve 


With the help of Equation 31 one can basically calculate the distribution function 
f for any time t > 0 if the initial distribution f(r,t = 0) is given. For any given 
initial distribution one has to use numerical methods, e.g. RUNGE-KUTTA or 
ADAMS-STORMER. Furthermore, one can ask for an analytical solution provided 
that the initial distribution is given by an analytical function. As a special case, 


10 


6 Equations for systems in which the initial distribution function is given by a 
narrow Gaussian curve 


assume that f(r,t = 0) is given by a narrow GAUSSian curve: 


(= SGserae (36) 


In Equation 36 Co is a constant, 79 is the mean particle radius at time t = 0 and 
Egro is the scatter of the initial distribution, where ¢) must be < 1 so that the 
GAUSSian curve is narrow. 
If one integrates the distribution function from r = 0 to r = oo, one gets according 
to Equation 15 the initially existing total number Zp of particles per unit volume 
of the dispersion. With 7 

r—To 


V2-€0Fo 


as an integration variable one gets for ¢5 < 1: 


U= 


i Pet = 0) dr = V2€070Co é i e ” du 
a ae 
+00 j 
= V2€ 9% oCo & Hl Ee iS) 2TEQTOC = Zo F (37) 


This gives the value for the constant Co in Equation 36: 


Zo 
V 20 - Eoro 


We then examine how the distribution function f(r,t) changes for t > 0. Initially, 
only the behavior within sufficiently short times is considered, during which the 
distribution function gradually broadens, but still remains relatively narrow (cf. 
Figure 1). Then one can assume that f(r,t) also for t > 0 is represented by a 
Gaussian curve and can be approximated with the spread e(t)r9 > oo, where 
e(t) <1. Accordingly it is set: 


C= (38) 


f(r,t) =C()- oo} (a8) (39) 


Since for e(t) < 1 the proportion of particles with molecular dimensions is 
vanishingly small, initially no particles disappear through complete dissolution. 
Therefore, the number of particles per unit volume remains the same Zp). With an 
equation analogous to Equation 37 we get: 


Zo 


NO Tn 


(40) 
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6 Equations for systems in which the initial distribution function is given by a 
narrow Gaussian curve 


| 


Figure 1: Distribution function after a given time t, provided that the distribution 
function for t = 0 is given by a narrow GAUSSian curve 


From Equation 39 and Equation 40 it follows: 
Zo _ 


oe eee eee 
f (r,t) ae ro) (41) 


In the case of a narrow GAUSS curve one can set r* = 79 in Equation 31, and 


also r = 79 , insofar as r appears as a factor. This gives us a simplified form of 
Equation 31: 


Of (r,t) 200eV kG 0 


Ot vRATH,+D or [(r — Fo) - f(r, 8] - (42) 


Substituting Equation 41 into Equation 42 gives the differential equation for <(t): 


de(t)  20aqV2 k@ 
dt vwrATkryot+D 


e(t). (43) 


Since r no longer appears in Equation 43, the heuristic approach in Equation 39 is 
justified. 
Taking into account the initial condition ¢(t = 0) = é9, integrating Equation 43 
results in: 

e(t) =ep-er (44) 
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6 Equations for systems in which the initial distribution function is given by a 
narrow Gaussian curve 


wherein for abbreviation 7 was set to 


—2 wins 
V5 AT kro +D 
= ; 45 
’ 20@V2 kD (45) 


The following limiting cases are to be noted: 
1. If kro >> J, i.e. for diffusion-controlled mass transfer it follows: 


vTrBART 


20 DoV,2 Be) 


TD 
According to Equation 46, the time constant Tp for diffusion-controlled mass 
transfer is proportional to the third power of the average particle radius 
79, and therefore grows disproportionately quickly with increasing average 
particle radius. A double mean radius corresponds to eight times the value of 
the time constant. Basically, the speed of dissolution-reprecipitation is always 
finite. In practice, however, the dissolution-reprecipitation is only noticeable 
when the time constant 7p is of the order of magnitude of the observation 
time t,p, or smaller. If 7p is larger than tops, then the size distribution can 
be considered practically stable for most purposes. Thus, from Equation 46 
one obtains as a condition for practical stability of a size distribution: 


F 205 
eat {222s tobs | (47) 


2. If kro < FY, i.e. for reaction-controlled mass transfer, it follows: 


a2 
VI5 AT 
ee eee 48 
‘i 2akcoV,2 (48) 


According to Equation 48 7p is only proportional to the square of the mean 
droplet radius. In view of the weaker dependence for Tg as a function of 7 
compared to the corresponding dependence for Tp according to Equation 46, 
a formula analogous to Equation 47 can be written down with only a lower 
degree of approximation. Despite this limitation, the condition for the 
practical stability of a size distribution for a given observation time is given 
analogously to Equation 47: 


Tos 


~N 


ee ‘Laby 


VAT a) 
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7 The behavior of a polydisperse system over long periods of time with 
diffusion-controlled mass transfer 


7 The behavior of a polydisperse system over long 
periods of time with diffusion-controlled mass 
transfer 


If the distribution function is initially given by a narrow GAUSS curve, the distri- 
bution curve will be broader with increasing time and the distribution function 
also assumes finite values for r = 0, so that particles according to Equation 16 
disappear. 

From Equation 30 one obtains for the time change in the radius of a single particle 
for diffusion-controlled mass transfer (kr >> Y) for r< r*: 


20 DoV? 
=f SS 50 
‘ vr? BAT ee) 
Inserting Equation 50 in Equation 16 yields: 
Z 2 ae pet 
og Se EDV, 0st) (51) 


dt r=0 vr? BT 


Since the right hand side of Equation 51 has a certain finite value under the 
conditions mentioned, it follows that f (r,t) increases proportionally to r? for r & 0. 
It is tentatively assumed that for longer observation times the distribution function 
f(r,t) can be represented as the product of a pure time function gp(t) and a 
function of the radius ratio Fe 
= By. 
¢ re (t) ( ) 
where r* is in turn a function of time. Taking into account the proportionality of 
f with r? for small values of r, let: 


f (r,t) = gno’hp(e), (53) 


where as normalization 


hp(g@=0)=1 (54) 


should apply. 
With regard to Equation 32 and Equation 53, the function hp(g) must satisfy the 
condition 


| hp(o)- 0 do = | hp(e) - 0° do. (55) 
0 0 
Substituting Equation 53 into Equation 17 gives by observing Equation 52: 


PM Tee) en nae, 68) 
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7 The behavior of a polydisperse system over long periods of time with 
diffusion-controlled mass transfer 


By inserting Equation 50, Equation 52 and Equation 53 into Equation 16 and 
eliminating gp(t) with the help of Equation 56, it follows: 


x 2 2 
dr - Oo 7 coVin (57) 
dé 38rqp [r*(H)] AT 
where for abbreviation: sz, 
YD = | hp(o)o" do. (58) 
0 


If the size distribution as a function of r is already given at time t = 0 by the 
general form according to Equation 53, then by integrating Equation 57 it follows 
that: 

t 


Tb) = = Oe es (59) 
Tp 
where 7 is a constant with the dimension of a time: 
, _ vplr(t=0R @T 
= : 60 
7D 20 DoV,2 ( ) 
Using Equation 56 and Equation 59 one gets from Equation 53: 
2, h 
POSE = o hole) - const. (61) 
3 


(+4) 
D 


Inserting Equation 61 into Equation 33 gives, taking into account Equation 52, 
Equation 59 and Equation 60 the differential equation: 


Ce | : hp (2) | e -3) Dew. (62) 


do YD 


Since the time t no longer appears explicitly in Equation 62, the tentative approach 
in Equation 53 is justified for quasi-stationary conditions. To calculate the function 
hp(e) from Equation 62, it is necessary to know the parameter yp, whose value is 
determined by Equation 55. 

If one multiplies Equation 62 by o® and performs integration by parts between 
o = 0 and 9g, one gets: 


[(£ -30+3) -ohp(0) +9: [f to(e)- o°do~ f ho(e)- ao] = 0. (63) 


In Equation 63 the difference of the integrals for @ — oo vanishes with respect to 
Equation 55. This results in the equivalent condition for determining the parameter 
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7 The behavior of a polydisperse system over long periods of time with 
diffusion-controlled mass transfer 


Yp instead of Equation 55: 


lim [(£ —30+ 3) hoi = 0: (64) 


The expression 


appearing in Equation 62 to Equation 64 has a minimum for @ = ,/yp with an 
associated function value 3 — 2: ,/yp. If one heuristically sets the latter expression 
equal to zero, one gets 7p = 9/4 and hereby from Equation 62: 


dIn(hp) eo -3 a i _ go-4 (65) 
do 40° —30+3 3+ 0 Cao 


Integrating and taking into account Equation 54 gives: 


7 


3 a) 3 3 = 0 
(esi) =(—B) (Gh) 
2 


fit (c > 5] ~0. (66b) 


With Equation 66a and Equation 66b we have Equation 64 for both 9 = oo and 
for the practical upper limit 9 = 3/2 fulfilled. If, on the other hand, one assumes 
yp < 9/4 or yp > 9/4, then Equation 64 cannot be fulfilled. For yp < 9/4, the 
function hp(@) calculated from Equation 62 becomes inversely proportional to 0° 
for large values of 0 and the left side of Equation 64 tends to a finite limit. For 
Yp > 9/4, hp(@) approaches infinity at a finite value of g and Equation 64 can also 
not be fulfilled. Thus, yp = 9/4 is the only value for yp that satisfies Equation 64. 
With yp = 9/4 one gets from Equation 60: 


dlised Ov [r*(t = OP BAT 
ee 80 D oV,2 


(67) 
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7 The behavior of a polydisperse system over long periods of time with 
diffusion-controlled mass transfer 


Inserting Equation 66a into Equation 61 yields: 


1+4 
D 
3 
for 9 = =a = BS (68a) 
i 3 
f(r, t) =0 for o= 7 (b) 2 3? (68b) 


where r*(t) is given by Equation 59. Equation 68a and Equation 68b are identical 
to the result of a calculation by LIFSHITZ and SLYOZOV, |9, 10], but their derivation 
however does not appear to be completely satisfactory, since there is no correspond- 
ing relationship like Equation 55 as a basis for determining the parameter yp. For 
a discussion of further details of Equation 68a and Equation 68b one may refer to 
the works of LIFSHITZ and SLYOZzov, [9, 10]. 

From Equation 68a one obtains for the value of o, at which f or In(f) has a maxi- 
mum as a function of 0, the equation 160?+270—54 = 0 with the solution 9 = 1.135. 


f(r) 


0 0,5 1,0 1,5 


Figure 2: Quasi-stationary distribution curve for diffusion-controlled mass transfer 
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7 The behavior of a polydisperse system over long periods of time with 
diffusion-controlled mass transfer 


Figure 2 shows the course of the distribution function obtained in this way for 
diffusion-controlled mass transfer. 
For t > Tp one obtains from Equation 59 and Equation 67 with 7 = r* according 
to Equation 64 the approximate equation: 


ee 3/80 D coV,2 -t 
PO == 9-7 = age (69) 


if kr > D and t > 7p or F(t) > F(t = 0). 

Equation 59 assumes that the quasi-stationary distribution according to Equation 53 
and Equation 66a as well as Equation 66b already existed for t = 0. However, if the 
distribution function for t = 0 is given by a moderately broad GAUSS curve, then 
the initial distribution transforms, according to the formulas of section 6, within 
a few time intervals Tp = ; - Ty into the quasi-stationary distribution existing for 
longer test times, so that Equation 69 can also be used in this case. 

Furthermore, the particle number Z per unit volume of the dispersion can be 
calculated as a function of time. If according to Equation 53 the distribution 
function for t 2 0 is the product of a function of time and a function of 9 = 
CLO SELF then the average particle volume is proportional to 7°. Since the time- 
idevendent volume fraction y of the dispersed particles is equal to the product 
of Z and the average particle volume, Z is inversely proportional to 7°. So using 
Equation 59 it follows: 


Z(t = 0) 


Z(t) = esa, if ke > . (70) 


It is noteworthy that the process of dissolution-reprecipitation assumed above as 
the cause of the decrease in the number of particles results in the same formal time 
law that SMOLUCHOWSKI, [{11, 12], obtained for the reduction of the number of 
particles by aggregation of colliding particles in the so-called rapid coagulation. In 
the formula of SMOLUCHOWSKI, which corresponds to the time law of a second- 
order reaction, the time constant is inversely proportional to Z(t = 0), while the 
constant TH in Equation 70 according to Equation 67 is proportional to [r*(t = 0)}° 
but independent of Z(t = 0). 

Relationships analogous to Equation 69 and Equation 70 can already be found in 
the work of TODES, [15], but without calculating the quasi-stationary distribution 
function in Equation 68a and Equation 68b. 
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8 The behavior of a polydisperse system over long 
periods of time with reaction-controlled mass 
transfer 


For reaction-controlled mass transfer, the calculation must be carried out in the 
same way as in section 7. 

From Equation 30 for the time change of a single particle one obtains for reaction- 
controlled mass transfer (kr < Y) for r<r*: 


dr QakepV2 
ge ELL vie 
dt t vr kT a) 
Insertion of Equation 71 into Equation 16 yields: 
dZ .. [2ckeoV2 f(r, t) 
eae a : Z 
Ti = vr RT we) 


Since the right-hand side of Equation 72 is finite under the conditions assumed here, 
it follows that f(r,t) for r = 0 increases proportionally with r or @. Analogous to 
Equation 53, but here with the factor @ instead of 0”, let: 


f(r,t) = gr(t)ohr(o), (73) 
where again 9 = rH according to Equation 52 and in the sense of normalization 
hy(o =0)=1. (74) 


With regard to Equation 34, the function hg(e) must satisfy the condition 


i, hr(o) - 0° do = if hp(o) - 0° do. (75) 
0 0 
Inserting Equation 73 into Equation 17 gives, taking Equation 52 into account: 


dn [r*(t)]" + fo” Pr(o) - ode 


By substituting Equation 71 and Equation 72 into Equation 16 using Equation 52 
and eliminating gr(t) with the help of Equation 76 it follows: 


gr(t) (76) 


dr* ———_ 2okeo V2 
dt — 3v7R [r*(t)] AT 


(77) 
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with e 
m= f hr(e)- odo. (78) 
0 


If the size distribution as a function of r is already given by Equation 73 at time 
t = 0, then by integrating Equation 77 it follows: 


r*(t) =r*(t =0)- ee (79) 


where Tp is a constant with the dimension of a time: 


, _ Brym [ett = PAT 
a Aa keoV,2 CO 


Using Equation 76 and Equation 79 it follows from Equation 73: 


_ehr(o) - const. (81) 
+4) 
R 


Inserting Equation 81 into Equation 35 gives, taking into account Equation 52, 
Equation 79 and Equation 80 the differential equation: 


(2 — 604 6) , dinlo) | (7 -6) -hy(Q) = 0. (82) 


YR do YR 


f(r, t) =i 


Since the time ¢t no longer appears explicitly in Equation 82, the approach in 
Equation 73 is justified for quasi-stationary conditions. To calculate the function 
hr(e) from the differential Equation 82, the parameter yg that occurs in it must 
first be determined, which is given by Equation 75. If Equation 82 is multiplied by 
o” and integrated by parts between 0 = 0 and a, it follows: 


(2& -60+8) - ghp(o) — 18- | tnio) edo [mete oP ae =0. (83) 


In Equation 83 the difference of the integrals for @ — oo vanishes with respect to 
Equation 75. This results in the following condition for determining the parameter 


YR: 


lim (72 —6o+ 6) hn) =0. (84) 


The expression 
2 


2 
oo 60456 
YR 
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appearing in Equation 82 to Equation 84 has a minimum for 9 = a with one 
associated function value 6 — 27 If one heuristically sets the latter expression 


equal to zero, one gets ya = ; and hereby from Equation 82: 


dIn(hr) 20-6 5 6 


_ = = ; 85 
do 302-60+6 2-0 (2-0) 85) 


Integration and taking into account Equation 74 gives: 


5 
halo S2)= (52) eH, (86a) 


he(e 2 2) =0. (86b) 


With Equation 86a and Equation 86b we have Equation 84 for both @ = oo and for 
the practical upper limit 0 =2 satisfied. If, on the other hand, one assumes yp < : 
or YR > :, then Equation 84 cannot be fulfilled. For yg > : the function hr(g) 
calculated from Equation 82 becomes for large values of 0 inversely proportional 
to o° and the left side of Equation 84 tends to a finite limit. For yp < , hr(o) 
approaches infinity for a finite value of 0 and Equation 84 cannot be fulfilled either. 
Therefore, only ya = ‘ comes into consideration. 

Under the conditions considered here, one gets for the mean radius 7 by using 
Equation 73 and Equation 86a and introducing 


as the integration variable: 


So fr)-rdr _ fo he(o)- do _ 8, (87) 
Io f(r) dr Jobr(e)-odo 9 


F= 


With yr = : one obtains from Equation 80 and Equation 87: 


,_virt(t=OP AT _ (2) weeeolat (88) 


7™R okcoV,2 8 okeoV,2 
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Insertion of Equation 86a into Equation 81 yields: 


(+4) = 
TR 
for p= . <2, (89a) 
r 
F(t) =0 for o=-— 22. (89b) 
r 


From Equation 89a one obtains for the value of @ at which f or In(f) has a maxi- 
mum as a function of 0 the equation 40? — 4 = 0 with the solution 0 = 1. 


0 0,5 1,0 1,5 2,0 
r3 £ 


r 9goT 


Figure 3: Quasi-stationary distribution curve for reaction-controlled mass transfer 


——_ 


Figure 3 shows the course of the distribution function obtained in this way for 
reaction-controlled mass transfer. 
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For t >> Tp the approximate equation is obtained from Equation 79, Equation 87 


and Equation 88: 
7 e akcoV2 -t 
le . m 90 
7D) VAT ey) 


kr < Dandt > 7, or rt) > rit = 0): 
Furthermore, for Z one obtains from Equation 15 by substituting Equation 73 and 
Equation 76 using Equation 52: 


7 39 ___, Jo hn(o) - edo 
4x -[re(t)!”  f> he(e)- odo | 


(91) 


According to Equation 78, the integral in the numerator is yr = .. For the integral 
in the denominator, using 


6 
Y= — 
2-0 
as the integration variable, one obtains: 
2 ioe) 
64 d 
| hee) - odo = —-— + 2°e° | ee SED 7, (92) 

0 9 3 u 


Inserting Equation 79, Equation 87 and Equation 92 into Equation 91 results for 
longer observation times in: 


gue # vRT \?* 
ae eae okcoV2 -t 


if kF > DY and t > 7; or F(t) > F(t = 0). 

As a special case of a reaction-controlled mass transfer, TODES and KHRUSCH- 
CHEW have calculated the redistillation of mercury between droplets of different 
sizes assuming the HERTZ-KNUDSEN evaporation equation and have compared 
their measurements with those of LIALIKOv, [7]. In the work of TODES and 
KHRUSHCHEV there are analogous relationships to Equation 69, Equation 70 and 
Equation 93, however, without the quasi-stationary distribution function having 
been calculated according to Equation 89a and Equation 89b 


(93) 
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In connection with experimental questions, the Equation 47, Equation 69 and 
Equation 70 for diffusion-controlled mass transfer as well as the Equation 49, 
Equation 90, Equation 91 and Equation 93 for reaction-controlled mass transfer 
are particularly important. The interfacial free enthalpy o appears in these equa- 
tions, the value of which is not readily available for crystals. The free enthalpy 
of the interface of ionic crystals with respect to the vacuum or saturated vapor 
is available from lattice theoretical calculations, {1]. Especially for ionic crystals 
with monovalent cations and monovalent anions, values of the order of magnitude 
100 erg cm~? are obtained. Significantly larger values are obtained for ionic crystals 
with higher-valued ions. However, the free enthalpy of the interface between an ionic 
crystal and an adjacent saturated aqueous solution, which is relevant in the present 
context, is significantly lower due to the reduction in the electrostatically induced 
parts of the energy because of the high dielectric constant of water. According 
to considerations by KAHLWEIT, [5], one can assume that the relevant value of 
o is between 20ergcm~? to 40 erg cm~? at room temperature, i.e. approximately 
the same order of magnitude as the interfacial tension between organic liquids 
and aqueous solutions. It should be expressly pointed out that the value of o 
for the interface between an ionic crystal and an aqueous solution depends on 
the composition of the latter and can be reduced in particular by surface-active 
additives. An excess of one of the two types of ions contained in the crystal can 
have the same effect. 

In Equation 47 and Equation 69 for diffusion-controlled mass transfer, only the 3rd 
root of o occurs. An estimated value a is therefore sufficient to indicate with the 
help of Equation 47 whether a size distribution is practically stable, or to calculate 
with the help of Equation 69 the average particle radius r after a given aging time. 
Here Equation 47 and Equation 69 can be used complementary, i.e. depending on 
the initial value 7(t = 0) the preconditions to use one or the other equation are in 
general fullfilled. From Equation 69 it follows, in accordance with known empirical 
facts, that the higher the solubility co, the larger the mean particle radius achieved 
after a given aging time. 

In order to facilitate the application of the above considerations, the mean radii 

r(t) expected according to Equation 69 for different solubilities co and observation 
times t are shown in Figure 4 using the following representative values: vy = 1, 0 = 
30ergcm~?, V;, = 30cm? mol’, J = 1x 10-> cm? s~!, T = 298° K (corresponding 
to 25°C). 
Diffusion-controlled mass transfer is only a possible limiting case. Existing observa- 
tions indicate that in many systems the rate of mass transfer is reaction-controlled. 
First of all, it generally follows from Equation 30 that the time change 7 of the 
radius of a particle has the greatest value when kr >> Y, i.e. when there is nearly 
an equilibrium at the phase boundary. Accordingly, the temporal changes of 7 and 
Z also have maximum values when kr > Y. 
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7 al 


10 


2 


10 10 10 10 10° 10° 
c, (Moles per litre} —~ 


Figure 4: Mean particle radius 7 after a given time t as a function of solubility 
co for diffusion-controlled mass transfer according to Equation 69 for 
r(t) > 7r(t =0) 


Now, in general, the value of & is not available. Without knowledge of the ratio a 


Equation 69 or Figure 4 only gives an upper limit of the achievable mean particle 
radius within a given observation time t . Accordingly, Equation 70 gives a lower 
limit for the particle number Z. 

If in a system the experimental value of r(t) is significantly smaller than the value 
for diffusion-controlled mass transfer calculated from Equation 69, then one can 
conclude that the mass transfer is essentially reaction-controlled. 

As an additional criterion for distinguishing between diffusion- and reaction- 
controlled mass transfer, one can use the shape of the distribution curve f(r, t) for 
small values, since with diffusion-controlled mass transfer f is proportional to r?, 
but with reaction-controlled mass transfer it is increasing directly proportional to 
r. This criterion is independent of an estimate for a, but requires a sufficiently 
precise experimental determination of the distribution function, especially for small 
values of r, which is by no means easy. 

However, Equation 49 and Equation 90 cannot easily be used for predictions, since 
the relevant quantity appears to be the rate constant k of the mass transfer between 
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the dispersed phase and the surrounding solution. In this case one can calculate k 
from observed values 7(t) using Equation 90: 


81 FOPAT 
ke oeV2E ie Teh Sr ESO): (94) 


Experimentally, Z is easier and more precise to determine than 7. Therefore, it is 
more convenient to calculate k using Equation 93: 


vAT ~ : ; 
k= (aoe). if 20 «2e=0). (95) 


The use of Equation 95 is particularly suitable for relatively poorly soluble sub- 
stances where determining the rate of growth or dissolution using the usual macro- 
scopic methods is difficult. Investigations into the reduction of k by specifically 
adsorbed components in the solution (higher-value foreign ions, organic dyes) are 
likely to be of particular interest, even if statements using Equation 95 are not 
immediately clear, since on the right hand side of Equation 95 an estimated value 
of o should be used, which is also influenced by solution additives, but probably 
generally to a lesser extent than k. 

It is outside the scope of this work to discuss the test results on individual systems. 
However, one fact should be pointed out. It is known that a precipitate of AgCl 
ages much more quickly than a precipitate of BaSO,. The solubilities and diffusion 
coefficients for both substances are approximately the same. According to the 
double charge of the ions in BaSOu,, a slightly larger value of o is to be expected for 
this substance than for AgCl. If the mass transfer for both systems were controlled 
by diffusion, the mean particle radius 7 that can be achieved after a given time t 
should be at least as large for BaSO, as for AgCl. The opposite finding is explained 
by the fact that the mass transfer for BaSOy, is reaction-controlled, as was shown 
by KAHLWEIT, [4], by tracking the temporal change in the electrical conductivity 
of supersaturated BaSO, solutions with particle sizes in the range of 0.1 1m to 
1m, while under similar conditions a diffusion-controlled mass transfer was found 
for AgCl. 

In the above considerations it is assumed that the solubility according to the 
GIBBS- THOMSON equation and the growth rate are unique functions of the particle 
size. These assumptions are not always fulfilled. A larger crystal may contain 
more structural defects and therefore may have a higher chemical potential than a 
smaller crystal. This is the only way it is possible for new nuclei to form on the 
base of a large crystal during the cathodic deposition of silver from silver nitrate 
solution, as KOHLSCHUTTER and TORRICELLI, [6], have observed. Furthermore, it 
is known that the growth rate of a crystal can depend significantly on the number 
of dislocations, so that different values of k shall be used for different crystal 
individuals. It is therefore possible that deviations from the distribution function 
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shown in Figure 3 occur, especially in reaction-controlled mass transfer. 

For crystals that are distributed in a second solid phase as a matrix, i.e. especially 
in alloys, another factor to be taken into account is the occurrence of local stresses 
that occur during the growth of the dispersed particles and which are compensated 
only within a certain delay through diffusion or plastic flow of the atoms of the 
surrounding matrix. LIFSHITZ and SLYOZOV, [9, 10], have pointed this out and 
announced a more detailed investigation. 


I would like to thank Mr. W. SENDLER for his help in working out the solutions 
to the differential equations Equation 62 and Equation 82, in particular in deter- 
mining the parameters that occur in these equations. 
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